Analytic perturbation theory versus expansion in the Gross-Neveu model 
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The expansion can be successfully used to calculate the Green functions of the two- 

dimensional 0{2N) Gross - Neveu model. In parallel, the methods of analytic perturbation theory 
are also applied. Comparing the results of these two calculations at leading order, we report on the 
surprising agreement between them. 
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An analytic model for the QCD running coupling 
as{Q'^) has been proposed in [ij. Since then, this concept 
evolved to become the method known as analytic pertur- 
bation theory (APT) [2]. The model is essentially based 
on the old idea of combining the renormalization group 
(RG) summation of the perturbative series of ultra-violet 
(UV) logarithms with analyticity in the variable. This 
idea has been used previously in the context of the QED 
ghost-pole problem Q. Certain reservations are in or- 
der when applying this idea tojthe QCD case. Some of 
them have been pointed out in [1]. However, this discus- 
sion did not address the pertinent aspect of QCD which 
relates to chiral symmetry breaking (ChSB) in the infra- 
red (IR) region [4]. It is commonly accepted that this 
nonperturbative phenomenon takes place in QCD. 

Recently attempts have been undertaken to apply a 
second order Bethe - Salpetcr formalism to the evalu- 
ation of the meson spectrum in QCD [5]. The idea of 
broken chiral symmetry was taken into account through 
an approximate solution of the appropriate Dyson - 
Schwinger equations. It has been shown that the me- 
son spectrum is well described if the analytic running 
couphng aApriQ^) is used in the region < 1 GeV^ , 
i.e. it was assumed there that except for analyticity, the 
IR behaviour of this coupling does not require any ad- 
ditional nonperturbative corrections, caused for instance 
by the dynamical ChSB effect in QCD. 

The analyticity of the running coupling function is a 
general requirement applied to any field theory, including 
the theories with dynamical ChSB. The latter however 
can possess an interesting property: there are models [6^1 
where the analyticity of the running coupling is auto- 
matically fulfilled in the true ground state. On the other 
hand, the same coupling, being calculated in the false 
vacuum, possesses a Landau singularity in the IR region 
of Q^. At large momenta both calculations coincide. 

Let us imagine for a moment that one knows only the 
result of the false vacuum calculation. It is natural in this 
case to use the APT method. The corresponding analyt- 
icization procedure, based on the spectral representation 
of the Kallen - Lehmann type, gives rise to a nonper- 
turbative contribution, which "cures" the IR ghost-pole 
problem. This improvement in the low-energy part of the 



theory should nevertheless be consistent with the idea of 
spontaneously broken symmetry: improved calculations 
based on the false vacuum should not deviate much from 
the result obtained from the true ground state. If devi- 
ations turn out to be big, one should assume that ad- 
ditional nonperturbative contributions must still be con- 
sidered. On the contrary, if deviations are small the APT 
result can be a good approximation even for theories with 
dynamical ChSB. 

Certainly, it would be most tempting to use the APT 
result without introducing additional corrections due to 
ChSB, i.e., supposing that the nonperturbative correc- 
tion inherent in analyticity already emulates a major part 
of the ChSB effect, at least up to some low-energy scale. 
The attempt to extract the IR QCD coupling from the 
meson spectrum and to cornpare it with the APT result 
has been made recently in Q . 

Here we show that there is a sensible way to discuss this 
point also on pure theoretical grounds by applying the 
APT method to a theory which possesses crucial prop- 
erties of QCD and yet is solvable both in the UV and 
IR regions. To put it more precisely, following the sol- 
uble two-dimensional 0(2Ai")-symmetric Gross - Neveu 
(GN) model @, which is known to be a renormalizable 
and asymptotically free quantum field theory, we argue 
in this letter that there is actually a remarkable agree- 
ment between the running coupling obtained by the APT 
method and, independently, by the 1/iV nonperturbative 
expansion. 

One should elucidate two points here. First, we com- 
pare only the leading order results, believing that higher 
orders can be interesting if the first step happens to be 
successful. Second, to be precise one can speak only 
about accurate agreement in the region of Euclidean mo- 
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above the perturbative Landau sin- 



gularity. It would be too naive to expect that the sim- 
ple analyticization procedure can perfectly describe all 
details of the nonperturbative dynamics. Nevertheless, 
our calculations clearly show that this procedure is com- 
patible with the fact that fermions acquire their masses 
dynamically. The same conclusion has been made in Q 
for QCD. 

The effective coupling g^{Q^) of the four-fermion in- 
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teractions in the GN model, obtained after an RG re- 
summation of the leading UV logs in the naive vacuum, 
i.e., for massless fermions, is described by the formula 
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(1) 



The l.h.s. is fixed as iV ^ cx), this explains the presence 
of the factor N. Here Pq = l/27r, A is an RG invari- 
ant scale, A = /^exp[— l/(2/3og^)], and is the couphng 
constant renormalized at a scale /i (a momentum cut-off 
scheme is used; dimensional regularization gives essen- 
tially the same result). 

After imposing APT analyticity one obtains an ex- 
pression similar to that found for the QCD coupling 



''APT 



(Q^) (see Ref. [ij) at the one loop order in the 
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The last term eliminates the unphysical Landau singu- 
larity at — A2 from the perturbative RG improved 
result ll]). This singularity results from the assumption 
that fermions are massless in the vacuum. We do this on 
purpose. Our aim at this stage is to obtain the invariant 
couphng g^{Q^) by the APT method without detailed 
description of the phase state of the GN model. In some 
respect this resembles the QCD case where one obtains 
the analogues of eqs. ^ and ^ neglecting details re- 
lated to the ChSB effect. As a possible alternative we 
mention here a model for the QCD invariant coupling, 
which is based on the hypothesis of finite gluon and quark 
masses Q. 

At present the main properties of the GN model are 
well studied. In particular, it is known that the 
expansion is a rather firm method to find that fermions 
are massive in the true vacuum Indeed, as has been 
stressed first by Witten ^1^, the GN model possesses 
a low-temperature phase of the Berezinski - Kosterlitz - 
Thouless type In this phase, the fermions have mass 
M ^ 0, just as if there were true spontaneous chiral 
symmetry breaking. It is fortunate (although perhaps 
just a coincidence) that the model exhibits in this way a 
feature that would be present in higher dimensions. Here 
we would like to stress only two important facts in favour 
of the 1/N expansion: (a) it yields the correct spectrum 
for the model J^] ; (b) the exact S-matrix obtained in [l^ 
can be expanded in powers of 1/iV correctly reproducing 
the known terms of the 1 /N series (see also Ref. [H). It 
shows that the 1/A^ expansion is a well convergent series. 

Thus, the expansion provides a solid basis for the 
calculation of the fermion mass M, which results in M = 
A at leading order (both the exact and results for 
the ratio M/A can be found in Ref. |l3)- To obtain 
the RG invariant running coupling cj^jQ'^) one can use 
the mass dependent RG method . In the true ground 
state {M ^ 0) one has % H 
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FIG. 1: The effective coupling of four-fermion interactions 
NPog'^^^'' as a function of the ratio Q^/A^. The curves repre- 
sent the functions given by eq. ^ (short-dashed curve), eq. 
Q (solid curve) and the RG resummed perturbative result 
eq. IT} (long dashes). Curves ([2]), (O intercept the ordinate 
at 1 and 0.5 respectively. The vertical dashed line corresponds 
to the position of the Landau singularity at /A^ = 1. 



where X = ^/l + 4M^/Q^ . Both expressions, namely, 
eqs. © and ([3]), have the same UV asymptotics. 

Note also that the fermion mass can be simply found 
by considering the one-loop effective potential V{(j)c) of 
the GN model as a function of a classical scalar field 4>c 
which represents the fcrmion-antifermion bound state 
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This function has two stationary points, (j)c = and 
4>c = A, and a positive second derivative, V"{A) = N/tt 
(minimum), for the second case alone. These two so- 
lutions correspond to the massless (false vacuum) and 
massive (true vacuum) fermion states. 

The functions given by eqs. (H]), and ^ are pre- 
sented in fig. [1] Our first observation is that the two 
functions Ng^p^Q^) and Ng'^j'^iQ'^) are quite close 
to each other in the interval 1 < Q2/A2 < oo. In this 
domain the analytic correction agrees perfectly well with 
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FIG. 2: The relative deviation of two quantities g^^r^iQ'^) 
and gy^^\Q^), described by eq. ([5]), is shown as a function 
of the ratio Q^/A^. 

the well-established fact that the fermions have a nonzero 
mass in the GN model. There is a sizeable difference 
in the results only at small values of Q^, i.e., in the 
IR region < < A^. Here, for instance, the value 
Ngl^'-^\0) = (2/3o)"^ is smaller than Ngll}^{0) by a fac- 
tor of two, and gl/^^\o) — 5apt(^^)- ^^^^ region, the 
results of the analytic model should be taken with care. 
This is our second observation. 

To compare the expressions obtained for the strength 
coupling, we consider the ratio 

-2 (1) _2 (1) „ 
_ 9 APT - 9l/N _ 
rL(X) — l2~(T) ' ^ ~ ^ ' ^ ' 

9l/N 

which describes numerically the relative deviation of the 
APT result from the prediction at point x. The 
function R{x), shown in fig. [2l decreases from i?(0) = 1 
to i?(ll.l) = —0.053, where it reaches its minimum. The 
line i? = is an asymptote for R(x) as x ^ oo. Inside 
the region 1 < x < oo the difference in values between 
the two approximations amounts to \R\ ~ 5% at most. 

The GN model with U{1) chiral symmetry has a mass- 
less pseudoscalar state n in the spectrum. The RG im- 
proved two-point Green function of this state, AjJQ^), 
obtained at leading order in the expansion Q, de- 
pends on the invariant coupling 

^NA(;HQ')=[l + ^^ Ngl/^\Q^). (6) 

The pseudoscalar propagator develops a pole at = 0, 
as it should be. On the other hand, its UV behaviour 
is described by formula (P) in accordace with general re- 
quirements of chiral symmetry. 

Now, suppose that the strength ^^/^'(Q^) in ([6]) is re- 
placed by the APT coupling 5apt(*3^)- ^^om our previ- 
ous consideration one might expect to see the strongest 
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FIG. 3: The 1/N approximation for the inverse propagator 
of massless pseudoscalar (full curve) in comparison with the 
APT result (dashed curve) as functions of Q^/A^. 



effect from such replacement in the IR domain of Q^. 
However, the physical IR singularity of eq. ^ a,t — 
dominates at small momenta, < < A^, causing a 
substantial suppression of this effect. This renders the 
agreement between the two approaches compared in fig. 
13] even better than between the quantities in fig. [T] (To 
avoid the singularity at = we plot the inverse func- 
tion of ([6]) in fig. [3]). Indeed, the two functions given by 
eqs. ^ and ([3]), which may differ up to a factor of 2 in 
the region < < A^ , show almost the same behaviour 
after being multiplied by the term (1 + 4A^/Q^). There- 
fore, although the low-energy behaviour of given by 
the running coupling 5apt('3^) ^^^^ accurate as com- 
pared with 5^/^^(Q^), the analytic approach is still suit- 
able for a reasonable description of the massless pseu- 
doscalar boson propagator in the IR region, provided the 
ChSB effects are correctly implanted. 

To conclude, we have compared here the leading order 
results of the standard expansion of the GN model 
in the Euclidean region of momenta > with the cor- 
responding results obtained by the methods of analytic 
perturbation theory. These two alternative approaches 
lead to different expressions for the effective RG invari- 
ant coupling constant g^ '■'^HQ^)- the first case the 
couphng g1l"^\Q'^) is essentially the result which con- 
tains important information about the dynamics of ChSB 
at low energies. In the second case the analytic coupling 
5apt(*3^) obtained without special regard to the ChSB 
effect, respecting only the basic principle of analyticity. 
The comparison shows that the nonperturbative contri- 
bution in ^aptIQ^) corrects the RG improved pertur- 
bative result (at least up to the region of the Landau 
singularity at = h?) in surprising accordance with 
the general requirements of dynamical ChSB in the IR 
region. This is an interesting and a priori unexpected 
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feature. We believe that our result is generic for renor- 
malizable theories with asymptotic freedom and dynam- 
ical chiral symmetry breaking. Therefore, with due care, 
it can also be applied to QCD, i.e. provided that the 
effect of dynamical ChSB is indeed the essential mech- 
anism governing the non-pcrturbative dynamics close to 
= 0. If other details are important we expect new 
contributions beyound the GN result to arise. 
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